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We  present  an  efficient  and  numerically  stable  method  to  calculate  time-dependent,  laser-induced 
temperature  distributions  in  solids  and  provide  a  detailed  description  of  the  computational  procedure 
and  its  implementation.  This  study  combines  the  two-dimensional  heat  equation  with  laser-induced 
heat  generation  and  temperature-dependent  luminescence.  The  time-dependent  optical  response  of  a 
system  is  obtained  numerically  by  the  Crank-Nicolson  method.  This  general  model  is  applied  to  the 
specific  case  of  optical  refrigeration  in  ytterbium  (Yb3+)  doped  fluorozirconate  glass  (ZB LAN).  The 
laser-induced  temperature  change  upon  optical  pumping  and  the  respective  transient  luminescence 
response  are  calculated  and  compared  to  experimental  data.  The  model  successfully  predicts  the 
zero-crossing  temperature,  the  net  quantum  efficiency,  and  the  functional  shape  of  the  transients.  We 
find  that  the  laser-cooling  transients  have  a  fast  and  a  slow  component  that  are  determined  by  the 
excited- state  lifetime  of  the  luminescent  ion  and  the  thermal  properties  of  the  bulk,  respectively.  The 
tools  presented  here  may  find  application  in  the  design  of  a  wide  range  of  optical  and  optoelectronic 
devices.  ©  2010  American  Institute  of  Physics,  [doi:  10.1063/1.3277009] 


I.  INTRODUCTION 

Thermal  processes  in  luminescent  materials  are  a  critical 
factor  determining  the  performance  of  solid-state  optical 
refrigerators,1  laser  gain  materials,2  display  and  lamp 
phosphors,3  solid-state  lighting  devices,4-6  luminescent 
biomarkers,7-9  and  fiber  lasers  and  amplifiers.  Light-induced 
heating  also  plays  a  key  role  in  laser  ablation,10  photothermal 
therapy,11-13  laser-induced  damage  of  tissue,  and  thermal 
lensing.  While  the  fundamental  aspects  of  this  problem  are 
readily  comprehended  from  the  basic  heat  equation,  actual 
calculations  of  time-dependent  laser-induced  temperature 
gradients  are  more  challenging  and  often  require  advanced 
numerical  tools  and  considerable  computational  power.14,15 
The  concept  of  laser-induced  heating  in  general  has  been 
studied  extensively,16-18  but  no  comprehensive  model  of 
laser-induced  heat  diffusion  and  its  effect  on  luminescence  in 
bulk  solids  exists. 

This  paper  presents  an  efficient  and  numerically  stable 
method  to  calculate  time-dependent  laser-induced  tempera¬ 
ture  distributions  in  solids  and  provides  a  detailed  descrip¬ 
tion  of  the  computational  procedure  and  its  implementation. 
In  addition,  we  introduce  a  temperature-dependent  material 
property,  in  this  case  luminescence  from  a  rare-earth  ion,  and 
show  that  the  respective  transient  response  of  the  system 
upon  optical  pumping  is  predicted  correctly  by  the  model.  To 
our  knowledge,  this  study  is  the  first  to  combine  the  heat 
equation  with  temperature-dependent  luminescence  in  order 
to  predict  the  time-dependent  optical  response  of  a  system. 
The  tools  presented  here  may  find  application  in  the  design 
of  a  wide  range  of  optical  and  optoelectronic  devices. 
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As  an  example,  we  apply  the  method  to  optical  refrig¬ 
eration  in  the  Yb3+  doped  fluorozirconate  glass  ZBLAN 
(ZrF4-BaF2-LaF3-AlF3-NaF).  The  development  of  laser- 
cooling  materials  and  devices  has  made  significant  progress 
over  the  past  decade.1,19  The  focus  has  been  primarily  on  the 
study  of  a  variety  of  rare-earth  doped  materials  and  ways  to 
fabricate  them  in  the  exceedingly  high  purity  and  optical 
quality  required  for  laser-cooling  applications.  However,  a 
quantitative  description  of  the  dynamics  of  laser-induced 
cooling  is  still  needed  and  is  critical  for  the  characterization 
of  laser-cooling  materials  as  well  as  for  the  design  and  per¬ 
formance  optimization  of  actual  optical  cryocooler  devices. 
The  relaxation  of  excited  rare-earth  ions  in  solids  involves 
both  radiative  and  nonradiative  processes.  The  nonradiative 
processes  are  exothermic  in  most  rare-earth  doped  materials; 
that  is,  net  heat  is  deposited  into  the  host,  and  the  material 
heats  as  a  result  of  laser  excitation.  Both  radiative  and  non¬ 
radiative  processes  also  occur  in  laser-cooling  materials;  but 
here,  the  net  result  of  the  nonradiative  processes  is  endo¬ 
thermic,  and  the  material  cools  as  a  result  of  laser  excitation. 
Two-band  differential  luminescence  thermometry  (TBDLT), 
a  technique  that  measures  laser-induced  changes  in  the  lumi¬ 
nescence  spectrum  in  the  time  domain  to  infer  subtle  changes 
in  internal  sample  temperature,  builds  on  transient  laser- 
induced  temperature  changes,  and  serves  as  a  test  of  the 
model  developed  in  this  study.  We  show  that  the  model  cor¬ 
rectly  predicts  the  time  dependence  of  the  TBDLT  signal  and 
finds  that  the  transients  have  a  fast  and  a  slow  component 
that  are  determined  by  the  excited- state  lifetime  of  the  lumi¬ 
nescent  ion  and  the  thermal  properties  of  the  bulk,  respec¬ 
tively. 

In  Sec.  II,  the  two-dimensional  (2D)  heat  equation  is 
introduced,  laser-induced  internal  thermal  processes  in  a 
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two-level  system  are  incorporated,  and  a  formal  description 
of  TBDLT  is  presented.  Section  III  applies  the  model  to  laser 
cooling  in  ZB  LAN :  Yb3+  glass  and  compares  the  calculated 
performance  with  experimental  TBDLT  transients. 

II.  MODEL  OF  TIME-DEPENDENT  LASER-INDUCED 
HEATING  IN  SOLIDS 

Consider  a  rectangular  solid  sample  that  is  excited  by  a 
single-mode  laser  focused  into  the  center  of  the  bulk  mate¬ 
rial.  If  scattering  can  be  neglected,  laser-induced  radiative 
and  nonradiative  processes  will  occur  over  the  focused  range 
of  the  laser,  that  is,  over  2 Zr  in  the  longitudinal  direction, 
where  Zr=(^ztI\p  is  the  Rayleigh  range.  Temperature 
changes  in  the  longitudinal  direction  over  the  Rayleigh  range 
will  be  minimal,  and  the  three-dimensional  system  can  be 
approximated  by  a  2D  transverse  slab  with  thickness  2 Zr.  We 
will  first  present  the  well-known  2D  heat  equation,  describe 
the  numerical  method  for  solving  the  respective  differential 
equations  (Sec.  II  A),  and  subsequently  introduce  laser- 
induced  heating  (Sec.  II  B).  Section  III  C  introduces  a  formal 
description  of  TBDLT  that  will  be  used  as  an  example  and 
validation  of  the  model  developed  in  this  section. 

A.  The  2D  heat  equation  and  the  Crank-Nicolson 
method 

The  diffusion  of  heat  is  governed  by  the  heat  equation, 
which  follows  from  the  Fourier  law  and  conservation  of  en¬ 
ergy.  The  Fourier  law 

4>q  =  -K\T  (1) 

states  that  the  local  heat  flux  <j)q  is  proportional  to  the  tem¬ 
perature  gradient  VT.  Note  that  the  thermal  conductivity  k 
generally  varies  with  temperature  and  direction  in  aniso¬ 
tropic  materials,  in  which  case  k  becomes  a  tensor.  In  the 
following  we  shall  ignore  both  these  dependencies.  Assum¬ 
ing  that  no  work  is  being  performed,  the  change  in  internal 
energy  per  unit  volume  A  Q  is  proportional  to  the  change  in 
temperature  AT;  that  is, 

AQ  =  CppAT.  (2) 

We  now  assume  that  the  specific  heat  capacity  Cp  and 
the  mass  density  p  are  both  independent  of  temperature.  In 
the  absence  of  internal  heat  generation,  the  change  in  internal 
energy  must  be  accounted  for  entirely  by  heat  flux  across  the 
boundaries,  and  therefore  the  change  in  internal  heat  and  the 
heat  flow  across  the  boundary  must  be  equal.  This  yields  the 
heat  equation,  which  in  the  2D  case  is  given  by 

dT(x,y,t)  k 
dt  Cpp 

In  Eq.  (3),  T(x,y,t)  is  the  temperature  (in  K)  at  time  t  (in  s) 
and  location  (x,y)  (in  m),  k  is  the  thermal  conductivity  (in 
W  m-1  K-1),  p  is  the  mass  density  (in  kg  m-3),  and  Cp  is  the 
specific  heat  capacity  (in  J  kg-1  K-1).  The  heat  equation  is  a 
second-order  partial  differential  equation  (PDE),  specifically 
it  is  a  parabolic  PDE.  Note  that  this  equation  does  not  ac¬ 
count  for  internal  heat  generation,  which  is  added  as  a  source 
term  later  in  this  section  [Eq.  (7)]  when  the  equation  is 


^T(v,y,0  ^T(v,y,0 

d?  +  W2 
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FIG.  1 .  Definition  of  quantities  for  the  CN  method  in  a  2D  square  grid. 


solved  using  the  Crank-Nicolson  (CN)  formalism. 

The  heat  equation  can  only  be  solved  analytically  in  a 
few  cases  and  usually  must  be  evaluated  numerically,  espe¬ 
cially  in  2D  and  three-dimensional  problems.  Several  ex¬ 
plicit  and  implicit  numerical  methods  exist  for  this  purpose. 
Explicit  methods  calculate  the  state  of  the  system  Y(t+At)  at 
a  later  time  t+At  from  the  state  of  the  system  Y(t)  at  the 
current  time  t;  that  is,  Y(t+ At)=F[Y(t)].  Implicit  methods  on 
the  other  hand  solve  an  equation  that  contains  both  the  cur¬ 
rent  and  the  later  state  of  the  system;  that  is,  G[Y(t),Y(t 
+  Af)]  =  0.  Explicit  methods  are  easier  to  implement  than  im¬ 
plicit  methods;  however,  they  often  fail  because  the  PDEs 
tend  to  be  unstable  unless  At  is  chosen  to  be  extremely  small, 
which  makes  explicit  methods  slow  and  sensitive  to  round 
off  errors.  In  contrast,  implicit  methods  require  upfront  com¬ 
putation  that  is  usually  more  than  offset  by  their  advantages 
of  unconditional  stability  and  larger  time  steps. 

The  CN  scheme  is  an  implicit  method21  that  is  a  particu¬ 
larly  powerful  approach  for  numerically  solving  parabolic 
PDEs  such  as  the  heat  equation.  It  is  a  second-order  method 
that  is  implicit  in  time  and  that  is  numerically  stable.  In  the 
2D  case  on  a  uniform  Cartesian  coordinate  grid,  the  heat 
equation  in  the  CN  scheme  becomes21 


(1+2  (7 tlj  +  Ttlj  +  Ttjlx  +  7?;-\)  =  ( 1 


where  TT  is  the  temperature  at  time  step  n  and  at  grid  loca¬ 
tion  (/,./)•  /z  is  the  dimensionless  Courant-Friedrichs-Lewy 
(CFL)  number  for  the  2D  case  and  is  given  by22 


KAt 


(5) 


where  Ad  is  the  lattice  constant  of  the  Cartesian  square  grid 
that  comprises  AX  A  cells.  These  quantities  are  illustrated  in 
Fig.  1.  Equation  (4)  can  be  written  in  matrix  notation  as 

[fin+lir+1]  =  [^iri  (6) 


where  [T”+1]  and  [7”]  are  one-dimensional  matrices  of  length 
A2  containing  the  temperatures  at  locations  (ij)  for  time 
steps  n+ 1  and  n ,  respectively.  [jmn+l]  and  [/zn]  are  2D  matri¬ 
ces  of  size  A2X  A2  that  contain  the  1  ±  2  /z  and  ±jul/2  fac¬ 
tors  in  Eq.  (4)  as  well  as  the  boundary  conditions  (see  Sec. 
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II  A  2).  In  the  presence  of  internal  thermal  processes  (e.g., 
laser-induced  heating  or  cooling)  with  power  P,  AQ  =  PAt  of 
thermal  energy  will  be  deposited  into  the  material  during  the 
time  interval  At.  According  to  Eq.  (2),  this  will  raise  the 
temperature  by  AT-AQ/Cpp-PAt/Cpp.  The  temperature  at 
time  step  n  + 1  can  now  be  found  by  multiplying  Eq.  (6)  with 
the  inverse  of  matrix  [/zw+1],  and  adding  the  effect  of  this 
additional  heat  source  such  that 

[r+1]  =  +  [at*].  (7) 

\/mn+l]  and  \/mn]  are  constant  for  a  given  system  if  k  and  Cp 
in  Eq.  (5)  are  assumed  to  be  independent  of  temperature.  If 
that  is  a  good  approximation,  [/zw+1]_1[/zw]  can  be  computed 
up  front,  and  the  method  then  proceeds  efficiently  from  one 
time  step  to  the  next  by  the  simple  matrix  multiplication  of 
Eq.  (7)  (see  Fig.  1),  yielding  the  time-dependent  spatial  tem¬ 
perature  distribution.  Sections  II  A  1— II  A  3  describe  the 
structure  of  the  /z  matrices,  the  implementation  of  boundary 
and  initial  conditions,  and  the  choice  of  time  step  in  the 
numerical  evaluation  of  Eq.  (7). 

1.  Structure  of  the  pi-matrices 

The  /z- matrices  in  Eq.  (7)  are  sparse  and  have  a  band- 
diagonal  structure  with  nonzero  elements  only  on  the  diago¬ 
nal  and  on  two  diagonals  on  either  side.  The  bandedness  is  a 
result  of  the  heat  flow  being  local,  i.e.,  heat  only  flowing 
between  neighboring  cells.  A  sparse  matrix  with  a  band- 
diagonal  structure  is  computationally  easier  to  invert  than  a 
dense  matrix;  however,  we  have  not  taken  advantage  of  this 
property  in  the  numerical  implementation  used  in  this  study 
and  have  used  a  standard  matrix  inversion  routine. 

Let  us  denote  the  diagonal  and  off-diagonal  elements  of 
the  /z-matrices  in  Eq.  (4)  by  a0  and  aq,  respectively.  For  the 
[/zw+1]  and  [jmn]  matrices,  these  elements  then  become 

ao+1  =1  +  2,11,  a1+1  =  -fi/  2,  (8a) 

<2q  =  1  -  2/z,  a\  =  /z/2.  (8b) 

To  illustrate  the  structure  of  the  resulting  /z-matrices,  let  us 
consider  a  simple  2D  system  consisting  of  4  X  4  cells.  The 
right  hand  side  of  Eq.  (6)  can  be  explicitly  written  as  shown 
in  Fig.  2  using  the  definitions  in  Eq.  (8b).  The  matrix  struc¬ 
ture  of  the  left-hand  side  of  Eq.  (6)  is  identical  but  uses  the 
elements  defined  in  Eq.  (8a).  Note  that  the  /z-matrix  consists 
of  A  rows  of  A  rows;  i.e.,  it  has  one  row  for  every  cell  of  the 
AX  A  Cartesian  grid. 

2.  Boundary  conditions  and  initial  conditions 

One  can  choose  either  Dirichlet  boundary  conditions, 
which  hold  the  boundary  at  a  given  temperature,  or  Neu¬ 
mann  boundary  conditions,  which  hold  the  boundary  at  a 
given  heat  flow  rate  (e.g.,  imperfect  insulation).  Here  we 
choose  Dirichlet  boundary  conditions,  which  place  the  sys¬ 
tem  in  a  bath  with  infinite  heat  capacity  and  a  fixed  tempera¬ 
ture;  i.e.,  the  temperature  of  the  cells  around  the  perimeter  of 
the  AX  A  system  does  not  change.  Specifically,  these  are  the 
cells  in  the  top  row,  bottom  row,  left  column,  and  right  col¬ 
umn  of  the  AX  A  Cartesian  grid  (see  Fig.  1).  The  matrix 


7=1 


2 


3 


4 


1=1 

2 

3 

4 

r 

n  n 

w 

ao  a\ 

ai  ■  ■  ■ 

n  n  n 

n 

al  aQ  a\ 

■  a\  ■  ■ 

n  n  n 

n 

al  a\ 

.  .  ax  . 

n  n 

n 

a\  a0 

a”  ..  . 

ao  a  1  • 

a[  .  .  . 

n 

WWW 

n 

ai  • 

ai  a\ 

■  ai  ■  • 

n 

n  n  n 

n 

.  ax 

al  a\ 

.  .  a  j  . 

n 

n  n 

n 

ClX 

d\  <lQ 

Cli 

n  n 

ax  .  .  . 

a0  a\ 

.  .  . 

n 

WWW 

n 

al  .  . 

^1  rto  a\ 

d  l  ■  • 

WWW 

■  •  "i  • 

^1  a 0  ^71 

■  ■  ai  • 

n 

w  w 

W 

a  1  n0 

n 

w  w 

ai  .  .  . 

ao  a\ 

■  a\  ■  • 

ai  a\ 

W 

WWW 

•  •  "i  • 

rtj  rtg  Cl  j 

M 

n  n 

'  •  '  "l 

al  ao 

K  i 
Th 
T\i 

T\  i 

T^ 

T\> 

mN 

1  V 
T\ 2 
Th 
Th 

T\ 3 

n, 

Tl 4 
T% 
Tl i 


FIG.  2.  Structure  of  the  matrices  on  the  right  side  of  Eq.  (6)  for  the  example 
of  a  4  X  4  Cartesian  grid.  Zero  elements  are  represented  by  dots  for  clarity. 


elements  for  these  cells  have  to  be  modified  to  ensure  con¬ 
stant  temperature.  Specifically,  the  first  and  last  A  rows  of 
the  above  /z-matrix  (see  Fig.  2)  correspond  to  the  top  and 
bottom  rows  of  the  AX  A  Cartesian  grid,  respectively,  and 
the  respective  diagonal  elements  are  replaced  with  1  and  the 
respective  off-diagonal  elements  with  0.  In  each  of  the  other 
blocks  of  A  rows,  the  first  row  and  last  row  correspond  to  the 
left  and  right  columns  of  the  AX  A  Cartesian  grid,  respec¬ 
tively;  again,  the  respective  diagonal  elements  are  replaced 
with  1  and  the  respective  off-diagonal  elements  with  0  to 
ensure  constant  temperature.  With  these  boundary  condi¬ 
tions,  the  matrix  structure  of  the  right  side  of  Eq.  (6)  (see 
Fig.  2)  is  modified  to  the  structure  shown  in  Fig.  3. 

Initial  conditions  are  provided  in  the  temperature  vector 
I1]  -  for  all  grid  locations  (ij)  at  t=0.  This  initial  temperature 
distribution  can  be  arbitrary.  However,  it  is  important  to  note 
that  the  temperature  of  the  perimeter  cells  of  the  AX  A  Car¬ 
tesian  grid  will  be  held  constant  at  their  respective  I^j  value 
by  the  boundary  conditions  defined  above. 
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FIG.  3.  Structure  of  the  matrices  on  the  right  side  of  Eq.  (6)  for  the  example 
of  a  4  X  4  Cartesian  grid  and  using  the  Dirichlet  boundary  conditions  of  Sec. 
II A  2.  Zero  elements  are  represented  by  dots  for  clarity. 


Downloaded  14  Dec  201 1  to  64.106.63.195.  Redistribution  subject  to  AIP  license  or  copyright;  see  http://jap.aip.org/about/rights_and_permissions 


063108-4  Patterson  et  al. 


J.  Appl.  Phys.  107,  063108  (2010) 


3.  Choice  of  time  step 


The  numerical  stability  of  the  CN  method  does  not  de¬ 
pend  on  the  size  of  the  time  step  Af.21  However,  there  is  an 
upper  limit  for  A t.  Note  in  Eq.  (8)  that  the  diagonal  elements 
should  be  sufficiently  close  to  1 .  Rather  than  choosing  A t,  it 
is  preferable  to  fix  /z  and  calculate  A t  from  Eq.  (5);  that  is, 


At- 


/xpCJM]2 


(9) 


With  proper  choice  of  /z,  this  approach  guarantees  numerical 
stability  and  the  maximum  possible  At  for  the  given  system. 
The  range  of  /z= 0.1,..., 0.2  has  proven  to  be  a  practical 
choice. 


B.  2D  heat  equation  with  laser-induced  heating  or 
cooling 

In  order  to  calculate  [7™+1]  in  Eq.  (7),  we  must  first 
obtain  the  rate  of  internal  heat  generation  P.  Assume  a  two- 
level  system  that  is  optically  pumped  and  that  can  decay 
radiatively  and  nonradiatively.  The  fraction  s  of  the  laser 
excitation  energy  is  converted  to  heat  (nonradiative  relax¬ 
ation)  while  the  fraction  1  -  s  undergoes  radiative  relaxation 
as  either  stimulated  or  spontaneous  emission.  Note  that 
stimulated  emission  is,  by  definition,  resonant  with  the  laser 
and  does  not  generate  any  heat.  The  rate  of  internal  heat 
generation  P(x,y,t )  is  therefore  proportional  to  the  sponta¬ 
neous  emission  rate  according  to 

he 

P(x,y,t)  =  n2(t)ANaVs  — ,  (10) 

Kp 

where  n2{t)  is  the  upper  level  population,  A  is  the  spontane¬ 
ous  decay  rate,  and  Na  is  the  number  density  of  absorbers.  In 
the  CN  scheme,  the  excited  volume  V  in  Eq.  (10)  corre¬ 
sponds  to  the  volume  of  one  cell,  V=2zr(AcI )2,  where  2 zr 
=  2irwl/\p  is  the  confocal  parameter  for  a  Gaussian  beam.23 

Following  the  analysis  of  a  two-level  system  allowing 
for  saturation  by  van  Dijk,  the  upper  level  normalized 
population  is  given  by 

n2  (*)  = - -{1 -exp[-(gwfl  +  A)f]},  (11) 

gwa+A 

for  continuous  constant  pumping  and  the  initial  condition 
n2( 0)  =  0.  In  Eq.  (11)  we  have  made  use  of  the  fact  that  the 
Einstein  coefficients  for  absorption  ( Bn )  and  stimulated 
emission  (#21)  are  related  by  B2l/ Bu=gl/ g2  and  have  intro¬ 
duced  g=  1  +  {g\l  g?),  where  gx  and  g2  are  the  degeneracies  of 
the  lower  and  upper  levels,  respectively.  The  degeneracy  is 
completely  lifted  in  the  low  coordination  symmetries  of  the 
glasses  studied  here,  and  we  therefore  assume  g\=g2  for  our 
computations.  In  Eq.  (11),  the  absorption  rate  constant  wa  is 
given  by 

wa(\,T)  =  (12) 

he 

where  IP(x,y,t )  is  the  laser  irradiance  (in  W  m-2)  and 
c ra(X,T )  is  the  absorption  cross  section  (in  m-2).  Note  that  in 
steady  state  and  for  high  irradiance  (7P^oo)  we  have 
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FIG.  4.  Schematic  of  energy  levels  in  Yb3+-doped  ZBLAN  showing  the 
ideal  case  of  laser  cooling  occurring  between  the  two  multiplets  of  Yb3+. 
The  dopant  is  excited  by  the  pump  laser  from  the  top  of  the  ground  state 
multiplet  (2F7/2)  to  the  bottom  of  the  excited- state  multiplet  (2F5/2).  Vibra¬ 
tional  energy  is  absorbed  from  the  host  during  thermalization  in  both  mul¬ 
tiplets,  and,  as  a  result,  the  material  cools. 


and  n2^l/2;  i.e.,  at  most  half  of  the  absorbers  can 
be  excited.  This  properly  accounts  for  saturation  and  is  con¬ 
sistent  with  the  fact  that  sustained  inversion  cannot  be 
achieved  in  a  two-level  system  under  adiabatic  conditions.25 
The  irradiance  at  a  CN  cell  at  location  (v,y)  is  given  by 

E(x,y,t)  =  P0(t)j^A,  (13) 

(A  dY 

where  the  incident  laser  power  P0(t )  (in  watts)  is  chosen  to 
be  a  step  function  at  t=t0.  The  normalized  transverse  spatial 
laser  power  distribution  <fi(x,y)  in  Eq.  (13)  is  assumed  to  be 
Gaussian;  i.e., 

<t>{x,y)  =  exp(_  2r2lw10),  (14) 

7TW0 

where  r2=x2+y2.  Equation  (14)  places  the  center  of  the 
Gaussian  beam  at  the  origin.  To  ensure  energy  conservation, 
the  physical  size  of  the  CN  grid  has  to  be  chosen  such  that  it 
substantially  contains  the  Gaussian  beam  profile;  i.e., 

rNAd/2 

<Kr)dr~  1.  (15) 

do 


C.  Model  of  differential  luminescence  thermometry  in 
solid-state  optical  refrigerators 

By  far,  the  most  widely  studied  solid-state  optical  refrig¬ 
erator  system  is  Yb3+  doped  into  the  fluorozirconate  glass 
ZBLAN.1  The  schematic  in  Fig.  4  illustrates  the  concept  of 
solid-state  laser  cooling  using  Yb3+.  A  pump  laser  is  tuned  to 
a  wavelength  ( Xp )  that  is  longer  than  the  mean  luminescence 
wavelength  (Xy),  and  the  energy  difference  corresponds  to 
the  amount  of  heat  that  is  extracted  as  heat  from  the  solid  for 
each  excitation/emission  cycle.  The  respective  laser-cooling 
efficiency  of  this  ideal  case  is  given  by  ycooi=(Xp-Xf)/Xf. 
Note  that  the  thermal  population  of  the  crystal  field  levels  in 
the  excited  state  is  temperature  dependent,  causing  Xf  to  red 
shift  and  rjcoo[  to  decrease  as  temperature  decreases. 
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In  a  real  system,  the  excited  state  may  have  a  net  quan¬ 
tum  efficiency  y(T,Xp),  which  is  less  than  unity  because  of 
nonradiative  processes  occurring  as  a  result  of  interactions 
with  impurities  as  well  as  background  absorption  due  to  di¬ 
rect  absorption  of  pump  energy  by  the  impurities.  Note  that 
77  is  a  function  of  both  temperature  and  pump  wavelength  not 
only  due  to  its  dependence  on  the  resonant  absorption,  but 
also  due  to  the  spectral  overlap  between  the  rare-earth  ion 
emission  and  the  impurity  absorption.  Incorporating  these 
effects,  the  cooling  efficiency  thus  becomes 


Xf(T) 


(16) 


As  the  temperature  is  lowered,  Xf  red  shifts,  rjcoo[  decreases, 
and  there  is  a  temperature  at  which  laser-induced  heating  and 
laser-induced  cooling  are  exactly  balanced;  i.e.,  T7cool=0. 
This  temperature  is  referred  to  as  the  zero-crossing  tempera¬ 
ture  (rzcx),  and  it  can  be  used  as  a  relative  measure  of  the 
concentration  of  impurities  and  the  overall  quality  of  a  laser¬ 
cooling  material.  The  net  quantum  efficiency  at  Tzcx  can 
thus  be  calculated  from  y(TZCT,Xp)  =  Xf(TZCT)/Xp.  At  tem¬ 
peratures  above  Tzcx,  ?7cool  >  0  and  heat  is  extracted  from  the 
solid,  while  below  ^ZCT»  77cooi  <  0  and  there  is  net  heating  of 
the  solid.  The  corresponding  laser-induced  rate  of  cooling  or 
heating  is  then  given  by  Eq.  (10).  The  factor  s  in  Eq.  (10) 
corresponds  to  the  negative  cooling  efficiency  (i.e.,  heat  gen¬ 
eration)  of  Eq.  (16).  Thus,  in  the  context  of  optical  refrigera¬ 
tion  of  solids,  laser-induced  internal  cooling  is  simply  con¬ 
sidered  as  negative  heating,  s  =  -ycoo i- 

Laser-induced  temperature  changes  inside  a  solid  can  be 
measured  by  observing  the  subtle  changes  that  occur  in  the 
material’s  luminescence  spectrum  as  the  laser  is  turned  on. 
This  technique  is  known  as  differential  luminescence  ther¬ 
mometry  (DLT);  specifically,  we  have  developed  TBDLT  as 
a  sensitive,  noncontact  method  to  characterize  laser-induced 
temperature  changes  in  solids.20  TBDLT  infers  changes  in 
the  local  sample  temperature  from  changes  in  the  lumines¬ 
cence  spectrum  that  occur  during  the  laser-induced  cooling 
process,  and  it  provides  a  temporally  and  spatially  resolved 
temperature  measurement  that  allows  for  rapid  performance 
screening  of  laser-cooling  samples.  Here,  a  brief  description 
of  TBDLT  is  presented,  and  it  will  be  used  in  Sec.  Ill  to 
validate  the  model  developed  in  Secs.  II A  and  II  B. 

The  luminescence  spectrum  of  Yb3+  in  ZBLAN  glass  is 
at  a  wavelength  near  1  yam.  As  shown  in  Fig.  5,  it  consists  of 
two  wavelength  regions  for  which  the  luminescence  intensity 
increases  (A  and  C)  and  two  wavelength  regions  for  which 
the  luminescence  intensity  decreases  (B  and  D)  as  the  tem¬ 
perature  is  raised.  Regions  A  and  D  are  spectrally  fairly 
broad  and  can  be  easily  selected  by  commercial  bandpass 
interference  filters.  The  TBDLT  method  chooses  these  two 
bands,  A  and  D,  and  detects  changes  in  their  relative  inten¬ 
sity  to  obtain  a  measure  of  the  associated  internal  tempera¬ 
ture  change.  At  a  given  location  (*,y)  with  temperature  TXJ, 
the  TBDLT  signal  is  defined  as 


+A  -B  +C  -D 


Wavelength  [nm] 

FIG.  5.  2F5/2^2F7/2  luminescence  spectra  of  Yb3+  in  ZBLAN  glass  at  vari¬ 
ous  temperatures.  The  gray  areas  indicate  wavelength  regions  in  which  the 
luminescence  intensity  increases  (A  and  C)  and  regions  where  the  lumines¬ 
cence  intensity  decreases  (B  and  D)  as  the  temperature  is  raised.  Portions  of 
regions  A  and  D  can  be  selected  with  commercially  available  interference 
Liters. 


IA(Txy,t)  -  fD(Tx  y,t ) 


(17) 


The  luminescence  intensities  IA(Txy,t)  and  fD(Txy,t )  in  Eq. 
(17)  are  integrated  over  the  product  of  the  luminescence 
spectrum  I(X,T )  and  the  bandpass  filter  transmission  spec¬ 
trum  0(\);  i.e., 


Il(Tx,y)  =  J  I(KTxJOA(k)dk, 

r„(Tx,v)  =  j  Kk,  TXJ  eD(k)dk .  (18) 


The  experimental  implementation  of  the  TBDLT  method 
uses  a  gain-balanced  amplified  pair  of  photodetectors  for  the 
simultaneous  measurement  of  IA(T),  I*D(T ),  and  I*A{T) 
-7*(T),  and  it  requires  that  the  optical  powers  reaching  the 
two  detectors  are  balanced.  In  practice  this  is  achieved  by 
simply  detuning  the  optical  alignment  of  the  band  with  the 
larger  signal  to  match  the  signal  of  the  other  band.  Here,  let 
us  introduce  a  factor  £  by  which  we  can  scale  I*D{T)  such  that 
it  matches  IA(T)  at  the  bath  temperature  T0;  i.e., 

JI(X,T0)QA(X)dX 
*  fI(KT0)0DQs)dX' 


Finally ,  since  luminescence  is  collected  from  the  entire 
pumped  volume,  we  must  account  for  the  temperature  distri¬ 
bution  in  the  pumped  volume  not  being  uniform.  The  mea¬ 
sured  TBDLT  signal,  H (T,t),  is  therefore  obtained  by  inte¬ 
gration  of  g(Txy,t)  [Eq.  (17)]  over  the  transversal  plane. 
Note  that  the  luminescence  intensity  is  proportional  to  the 
excited-state  population  n2(x,y,t)  [Eq.  (11)],  and  £( Txy,t )  is 
thus  weighted  by  n2(x,y,t)  according  to 


S  (T,t)  = 


n2(x,y,  t)  £(  Txy,  t)dxdy , 


(20) 


where  T  becomes  the  average  temperature  in  the  optically 
excited  volume.  With  these  definitions,  we  obtain  E (T 
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<r0)<0,  H(r=r0)  =  0,  and  H(r>r0)>0.  That  is,  the  TB- 
DLT  signal  will  become  positive/negative  upon  laser-induced 
heating/cooling  of  a  system  that  was  initially  thermalized  at 
To¬ 
ni.  RESULTS  AND  DISCUSSION 

In  the  following,  the  model  developed  in  Sec.  II  is  ap¬ 
plied  to  optical  refrigeration  in  Yb3+-doped  ZBLAN  glass  as 
both  an  example  and  a  quantitative  validation.  In  Sec.  Ill  A 
we  will  first  present  the  relevant  material  parameters  of 
Yb3+-doped  ZBLAN  glass  that  are  required  for  the  calcula¬ 
tion  of  laser-induced  temperature  changes  (Sec.  Ill  B)  using 
the  model  developed  in  Sec.  II.  Section  III  B  presents  a  cal¬ 
culation  of  the  laser-induced  temperature  distribution  in 
ZBLAN :  Yb3+  as  a  function  of  time,  and  Sec.  Ill  C  shows  the 
transient  response  of  the  respective  TBDLT.  Finally,  a  dis¬ 
cussion  of  sample  size  effects  is  given  in  Sec.  Ill  D. 

Two  Yb3+-doped  fluorozirconate  glass  samples  were 
used  for  the  measurements  in  this  study.  Sample  I  was  a 
ZBLANI :  l%Yb3+ 

(ZrF4-BaF2-LaF3  -  A1F3  -NaF-InF3  -  YbF3)  glass  fabri¬ 
cated  in  our  laboratory  from  purified  precursor  materials. 
Sample  II  was  a  ZBLAN:2%Yb3+  glass  obtained  from  IPG 
Photonics  and  known  to  be  a  good  laser  cooler.  The  samples 
were  mounted  in  a  liquid-nitrogen  flow  cryostat  in  close  ther¬ 
mal  contact  with  the  cold  finger.  The  sample  temperature  was 
measured  with  a  miniature  temperature  sensor  mounted  di¬ 
rectly  on  the  sample  surface.  Optical  excitation  was  achieved 
by  a  single-pass  geometry  with  focusing  the  laser  into  the 
center  of  the  sample. 

A.  Spectroscopic  and  material  properties  of 
ZBLAN  :Yb3+  for  thermal  diffusion  modeling 

Several  spectroscopic  parameters  of  ZBLAN  :Yb3+  are 
needed  for  the  thermal  diffusion  calculations.  The 
temperature-dependent  absorption  cross  section  at  the  pump 
wavelength  [Eq.  (12)]  was  obtained  by  fitting  a  cubic  poly¬ 
nomial  to  the  measurements  of  the  absorption  cross  section 
at  Xp=  1020.6  nm  at  different  temperatures  [see  Fig.  6(a)]. 
Likewise,  the  temperature  dependencies  of  the  mean  lumi¬ 
nescence  wavelength  Xy  [in  Eq.  (16)],  as  well  as  the  lumines¬ 
cence  intensity  integrals  JI(X,T)6A(X)dX  and 
JI(X,T)0D(X)dX  [Eq.  (18)],  were  obtained  by  fitting  cubic 
polynomials  to  respective  experimental  data  obtained  at  dif¬ 
ferent  temperatures  [see  Figs.  6(b)  and  6(c)].  Table  I  summa¬ 
rizes  the  respective  polynomial  coefficients.  Table  II  summa¬ 
rizes  the  material  properties  typical  of  ZBLAN  glass  as  well 
as  the  laser  and  computational  parameters  specific  to  the 
present  experiments. 

B.  Laser-induced  temperature  changes  in  ZBLAN:  Yb3+ 

The  thermal  response  [Eq.  (7)]  and  the  resulting  lumi¬ 
nescence  response  [Eq.  (20)]  of  a  solid  under  Gaussian  beam 
irradiation  can  now  be  calculated.  As  an  example,  Fig.  7 
shows  four  calculations  of  the  time-dependent  transversal 
temperature  distribution  (at  T0  =  300  K)  for  four  hypothetical 
samples  of  ZBLAN :  Yb3+  having  different  net  quantum  effi- 


FIG.  6.  Temperature  dependence  of  spectroscopic  parameters  of 
ZBLAN:Yb3+:  (a)  absorption  cross  section  cra(\,T)  at  Xp=  1020.6  nm  (Ref. 
26),  (b)  mean  luminescence  wavelength  Xf  derived  from  luminescence  spec¬ 
tra  at  different  temperatures,  and  (c)  spectral  overlap  integrals 
f  I(X  ,T)  0A(X)dX  (circles)  and  fI(X,T)6D(X)dX  (squares)  [see  Eq.  (18)].  The 
solid  lines  in  each  of  the  plots  are  least-squares  fits  to  a  cubic  polynomial, 
and  the  respective  parameters  are  summarized  in  Table  I. 

ciencies  (77=1,  0.98,  0.97,  and  0.94).  Continuous  wave  (cw) 
laser  excitation  with  a  Gaussian  beam  profile  begins  as  a  step 
function  at  t=0  in  these  calculations.  The  lower  part  of  Fig. 
7  shows  the  corresponding  change  in  temperature  at  the  cen¬ 
ter  of  the  Gaussian  beam  over  the  course  of  5  s.  The  top  row 
in  Fig.  7  represents  the  case  of  an  ideal  ZBLAN  :Yb3+ 
sample  that  has  77= 1  and  exhibits  laser-induced  cooling.  For 
this  case,  the  center  of  the  pumped  area  cools  from  300  to 
299.587  K  during  the  first  20  ms  and  reaches  a  steady-state 
temperature  of  299.18  K  (at  t=5  s).  Samples  with  lower  net 
quantum  efficiency  show  less  laser-induced  cooling  or  even 
laser-induced  heating,  as  illustrated  for  the  77=  0.94  case  (Fig. 
7,  bottom  row).  Note  that  internal  heat  generation  is  essen¬ 
tially  zero  at  77=0.97.  This  is  expected  from  Eq.  (16):  77cool  is 
zero  for  77=Xy/Xp  =  0.9718;  i.e.,  300  K  is  the  Tzct  for 
ZBLAN:  Yb3+  with  77=0.9718  pumped  at  1020.6  nm.  Net 
laser-induced  heating  thus  occurs  for  77c  0.97. 

C.  Transient  response  of  ZBLAN  :Yb3+  luminescence 

The  TBDLT  transients  measured  for  samples  I  and  II  at 
various  bath  temperatures  are  shown  in  Fig.  8.  Laser-induced 
cooling  is  evident  by  a  decrease  in  the  TBDLT  signal  H (T,t). 
The  dependence  of  the  TBDLT  signal  as  a  function  of  time  is 
linear  in  the  double-logarithmic  representation  of  Fig.  8,  in¬ 
dicating  that  it  follows  a  simple  power  law  according  to 

S(r)«r#.  (21) 

The  slope  #  of  the  line  in  double-logarithmic  representation 
is  a  metric  for  the  laser-induced  temperature  change,  and  it 
can  be  used  as  a  figure  of  merit  for  the  laser-cooling  perfor¬ 
mance  of  the  material  at  a  given  bath  temperature.  Laser- 
induced  cooling  or  heating  is  therefore  present  if  #<0  or 
#>0,  respectively.  Figure  9  presents  the  #  values  obtained 
from  fits  of  Eq.  (21)  to  the  experimental  data  (Fig.  8),  and 
Tzc t  was  estimated  to  be  158  and  238  K  for  sample  II  (filled 
squares)  and  sample  I  (filled  circles),  respectively.  At  these 
temperatures,  the  mean  luminescence  wavelengths  are  found 
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TABLE  I.  Coefficients  obtained  from  least-squares  fits  of  the  function  y(T)=a0+alT+a2T1+a3T 3  to  the  experi¬ 
mental  data  shown  in  Fig.  6. 


Coefficient 

Absorption  cross  section,  cra 
(m2) 

Mean  luminescence  wavelength,  \f 
(nm) 

Band  A 
(a.u.) 

Band  D 
(a.u.) 

a0 

1.29  X  10_2S 

1011.33 

-2.58 

0.669 

ax 

-3.90  X  10"28 

-0.130 

0.030 

0.071 

a2 

3.09  XIO-30 

3.2  X  10"4 

-4.64  X10"5 

-2.91  X  lO"4 

a3 

-4.46  X  10“33 

-3.52  X  10“7 

4.09  X  10“8 

3.73  x  10“7 

to  be  995.9  and  993.96  nm,  respectively  [see  Fig.  6(b)  and 
Table  II].  With  a  fixed  laser  excitation  wavelength  of  1020.6 
nm,  the  net  quantum  efficiencies  y  at  Tzct  are  thus  calcu¬ 
lated  to  be  0.9758  (sample  II)  and  0.9739  (sample  I)  [Eq. 
(16)].  Figure  9  also  shows  #  values  calculated  from  the 
model  presented  in  Sec.  II  [Eq.  (20)].  The  calculations  as¬ 
sumed  the  above  net  quantum  efficiencies,  and  a  scale  factor 
was  applied  to  the  experimental  data  to  match  the  absolute 
value  of  #  at  300.25  (sample  I)  and  257.05  K  (sample  II). 
The  model  provides  a  good  quantitative  description  of  the 
experimental  data  and  thus  serves  as  a  useful  tool  to  predict 
the  laser-cooling  performance  of  other  materials. 

An  interesting  feature  is  observed  at  low  temperatures 
for  sample  II.  #  and  thus  laser-induced  heating  reach  a  maxi¬ 
mum  at  133  K.  This  maximum  is  the  result  of  two  counter¬ 
acting  effects.  On  the  one  hand,  the  laser-cooling  efficiency 
gradually  decreases  with  decreasing  temperature  [Eq.  (16)] 
causing  more  of  the  absorbed  power  to  be  converted  to  heat 
[Eq.  (10)].  On  the  other  hand,  the  absorption  coefficient  at 
the  pump  wavelength  decreases  rapidly  with  decreasing  tem¬ 
perature  [Fig.  6(a)].  Therefore,  T— >0  causes  cra^ 0  causing 
the  rate  of  internal  heat  generation  P— >0  and  thus  a 

trend  that  is  confirmed  by  the  experimental  data. 


TABLE  II.  Summary  of  key  parameters  characterizing  the  ZB  LAN :  Yb3+ 
glass,  the  laser  excitation,  and  the  CN  computation. 


Parameter 

Symbol 

Value 

Units 

Ref. 

Material  parameters 

Thermal  conductivity 

K 

0.77 

W  m-1  K"1 

26 

Density 

P 

4445 

kg  m-3 

26 

Specific  heat 

c„ 

670 

J  kg’1  K-1 

26 

Specific  heat  capacity 

Cp=Cvp 

2.9782  X106 

J  m“3  K-> 

Yb3+  ion  density  (1  mol  %) 

Na 

1.899  X  1026 

m-3 

27 

Radiative  relaxation  rate 

MV 

540 

s-1 

27 

Internal  net  quantum  efficiency  77 

Varied 

Laser  parameters 

Laser  power 

% 

3.75 

W 

Laser  wavelength 

1020.6 

nm 

Beam  waist 

w0 

10 

yum 

Laser  on  time 

% 

0 

s 

CN  parameters 

Square  lattice  constant 

A  d 

50 

/xm 

Square  grid  elements 

NXN 

71X71 

CFL  number 

ix 

0.2 

D.  Sample  size  and  associated  characteristic  time 
constants 

The  model  also  allows  the  study  of  the  dependence  of 
the  transient  response  on  the  size  of  the  sample.  This  is  il¬ 
lustrated  in  Fig.  10,  which  shows  calculated  TBDLT  tran¬ 
sients  over  the  course  of  60  s  for  ZBLAN:Yb3+  with  y 
=  0.9739  and  various  2D  sample  sizes.  Three  regimes  are 
found.  There  is  an  initial  fast  component  that  is  independent 
of  sample  size  and  that  has  a  time  constant  governed  by  the 
excited- state  lifetime  of  the  rare-earth  ion.  The  respective 
time  constant  is  characteristic  of  how  quickly  heat  is  re¬ 
moved  from  (cooling)  or  deposited  into  (heating)  the  sample 
in  the  small  excitation  volume  defined  by  the  focused  laser. 
This  time  constant  is  on  the  order  of  only  a  few  milliseconds 
in  ZBLAN :  Yb3+  and  was  not  resolved  in  the  measurements 
shown  in  Fig.  8.  The  initial  fast  component  is  followed  by  a 
slower  component,  the  duration  of  which  depends  on  the 
sample  size.  The  time  constant  of  this  component  is  gov¬ 
erned  by  the  heat  capacity  and  thermal  conductivity  of  the 


t  -  20  ms  Is  3s  5s 


0.01  0.1  1 
Time  [seconds] 


FIG.  7.  Spatial  (top)  and  temporal  (bottom)  representations  of  calculated 
laser-induced  heat  diffusion  in  ZBLAN:  l%Yb3+  glass  pumped  at  1020.6 
nm.  The  calculation  was  performed  using  the  2D  CN  method  in  a  3.55 
X  3.55  mm2  grid  (see  Sec.  II)  and  the  parameters  of  Tables  I  and  II.  The 
internal  net  quantum  efficiency  77  is  reduced  from  top  to  bottom  in  the 
figure.  The  gray  scale  covers  the  range  of  299.0  K  (black)  to  300.1  K  (white) 
with  the  bath  temperature  held  at  T0=300  K. 
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0.55  x  0.55  mm2 
1.05  x  1.05  mm2 

1.55  x  1.55  mm2 
2.05  x  2.05  mm2 

2.55  x  2.55  mm2 
3.05  x  3.05  mm2 

3.55  x  3.55  mm2 


FIG.  10.  TBDLT  transients  for  ZB  LAN:  Yb3+  calculated  from  Eq.  (20)  with 
77=0.9739  and  varying  the  sample  size.  The  spatial  resolution  was  held 
constant  (50  X  50  jmm2  grid  element  size),  and  the  sample  size  was  varied 
via  the  grid  dimension  N.  Material  parameter  values  from  Tables  I  and  II 
were  used. 


first  few  seconds  are  well  described  by  assuming  a  sample 
size  larger  than  3X3  mm2  (see  Fig.  10),  which  is  consistent 
with  the  transversal  dimension  of  the  actual  samples  used  in 
our  experiments. 


FIG.  8.  Experimental  TBDLT  transients  of  ZBLAN:Yb3+  samples  I  and  II 
(see  Sec.  Ill  A)  recorded  at  different  bath  temperatures,  T0.  The  solid  lines 
are  fits  to  Eq.  (21),  and  the  respective  if- values  are  shown  in  Fig.  9. 

material  and  extends  for  several  seconds  for  ZBLAN:Yb3+ 
samples  with  cross-sectional  areas  of  >5  mm2.  This  is  the 
time  regime  in  which  the  measurements  of  this  study  were 
carried  out  (Fig.  8)  and  from  which  the  TBDLT  parameter  # 
was  calculated  (Fig.  9).  Finally,  the  thermal  processes  in¬ 
duced  by  a  cw  laser,  combined  with  the  sample  surface  being 
held  at  the  constant  bath  temperature  T0,  produces  a  steady- 
state  temperature  distribution  (and  thus  a  constant  TBDLT 
signal)  after  some  longer  time.  The  time  period  for  this 
steady  state  to  develop  in  larger  ZBLAN  samples  is  on  the 
order  of  many  seconds  and  is  determined  by  the  total  amount 
of  heat  being  deposited  into  the  sample  (i.e.,  sample  size) 
and  the  thermal  properties  (thermal  conductivity  and  heat 
capacity).  Note  that  the  initial  fast  temporal  response  of  the 
0.55  X0.55  mm2  grid  shown  in  Fig.  10  is  typical  of  what 
would  be  expected  for  a  sample  with  a  small  cross-sectional 
area  such  as  a  bare  optical  fiber.  In  the  case  of  bulk 
ZBLAN  :Yb3+,  the  measured  TBDLT  transients  during  the 


Temperature  [K] 

FIG.  9.  TBDLT  parameter  it  [Eq.  (21)]  for  experimental  data  (filled  sym¬ 
bols)  and  calculated  data  (open  symbols)  for  samples  I  and  II  (see  Sec.  III). 
The  Tzcx  is  found  to  be  238  and  158  K  for  samples  I  and  II,  respectively. 


IV.  CONCLUSIONS 

We  have  presented  a  quantitative  model  that  (1)  de¬ 
scribes  the  time-dependent  laser-induced  temperature  distri¬ 
bution  in  a  solid  and  (2)  correlates  the  laser-induced  tem¬ 
perature  changes  with  changes  in  luminescence  properties. 
The  implicit  CN  scheme  used  for  the  computational  evalua¬ 
tion  of  the  heat  equation  was  found  to  be  numerically  stable 
and  efficient,  allowing  for  the  rapid  exploration  of  the  pa¬ 
rameter  space.  As  an  example,  the  model  was  validated  for 
TBDLT  in  ZBLAN  :Yb3+  optical  refrigerator  samples.  The 
laser-cooling  performance  as  a  function  of  temperature  was 
accurately  predicted  by  the  model,  allowing  the  ^ZCT  and  net 
quantum  efficiency  to  be  calculated.  Furthermore,  the  model 
revealed  the  presence  of  three  distinct  time  constants  that 
govern  the  luminescence  response  in  optical  refrigerators 
upon  laser  excitation.  The  tools  developed  in  this  study  are 
general  and  can  be  readily  applied  to  other  materials  and 
temperature-dependent  properties,  enabling  quantitative 
studies  of  light-matter  interactions  in  a  wide  range  of  mate¬ 
rials  and  devices. 
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